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Abstract 

We use (3+1) dimensional ideal hydrodynamics to describe the space-time evolution of strongly 
^ '. 

interacting matter created in Au+Au and Pb+Pb collisions. The model is applied for the domain 
■ of bombarding energies 1-160 AGeV which includes future NICA and FAIR experiments. Two 



(N 



X 



equations of state are used: the first one corresponding to resonance hadron gas and the second 



^ I one including the deconfinement phase transition. The initial state is represented by two Lorentz- 

00 ■ 

■ boosted nuclei. Dynamical trajectories of matter in the central box of the system are analyzed. 

a^ : 

■ They can be well represented by a fast shock-wave compression followed by a relatively slow 
I isentropic expansion. The parameters of collective fiows and hadronic spectra are calculated under 

assumption of the isochronous freeze-out. It is shown that the deconfinement phase transition 
leads to broadening of proton rapidity distributions, increase of elliptic fiows and formation of the 
directed antifiow in the central rapidity region. These effects are most pronounced at bombarding 
energies around 10 AGeV, when the system spends the longest time in the mixed phase. Prom 
the comparison with three-fiuid calculations we conclude that the transparency effects are not so 
important in central collisions at NICA-FAIR energies (below 30 AGeV). 

PACS numbers: 24.10.Nz, 25.75.Ld, 25.75.Nq 
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I. INTRODUCTION 



Relativistic hydrodynamics is a very powerful tool to study high-energy nuclear collisions. 
Especially useful is its sensitivity to the "equation of state" (EoS) of strongly interacting 
matter and, in particular, to its phase diagram. In fact, extracting this EoS the main 
goal of heavy-ion experiments. Theoretical studies of the EoS are still not very conclusive. 
First principle QCD calculations on the lattice jl, 2| give reliable results only at small baryon 
chemical potentials. In this case a cross-over transition from hadronic to quark-gluon degrees 
of freedom is predicted at temperatures T ~ 170 MeV. Some signatures of a baryon-free 
quark-gluon plasma with low viscosity have been already found ^ in RHIC experiments with 
cm. bombarding energies a/sJvjv = 60 — 200 GeV. These signatures may be even stronger at 
LHC energies {^/SNN ~ 5 TeV). On the other hand, many phenomenological models predict 
that a strong first order phase transition may occur in compressed baryon-rich matter j4-6|. 
Presumably, such matter is created at lower (SPS, AGS) energies. A more detailed data 
should be obtained in the low-energy runs at RHIC as well as in future NICA [S] and 
FAIR jsl experiments. 

A large amount of experimental data for nuclear collisions at AGS, SPS and RHIC energies 
have been successfully described by different versions of the hydrodynamic model. The first 
model of this kind has been proposed by Landau more that 50 years ago |lol]. Unless stated 
otherwise, below we are dealing with perfect fluids, i.e. we neglect possible dissipative terms, 
associated with viscosity, heat conductivity as well as chemical non-equilibrium effects. In 
other words, it is assumed that deviations from local equilibrium are small starting already 
from early stages of a nuclear collision. One can roughly divide the existing versions of 
ideal hydrodynamics into two classes. The first class includes the models which apply fluid 
dynamical simulations from the very beginning i.e. starting from cold equilibrium nuclei. 
The attractive feature of such an approach is that no additional parameters are needed 
to characterize the initial state of the reaction. The models of the second class introduce 
an excited and compressed initial state - a locally equilibrated "fireball". It is believed that 
such fireball is formed at an early non-equilibrium stage of the collision. The disadvantage 
of this approach is a large freedom in choosing geometrical and fluid-dynamical parameters 
of the initial state. Due to this reason the predictive power of second-class models is greatly 
reduced, especially when studying sensitivity of the results to the EoS. Up to now many 
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versions of the fireball-based hydro dynamic model were developed ranging from simplified 
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dimensional ones 
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- Il6| and (2+l)-dimensional models [13|, ll7H20l| to more sophisticated (3+1)- 



26] ■ Recent calculations with inclusion of dissipative terms 27|, |28| 



show that data at RHIC energies can be reproduced with rather low viscosity coefficients. 
Historically, early 3D models of relativistic nuclear collisions 
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34| used cold Lorentz- 
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40j. The most 



contracted nuclei in the initial state. It is believed that such models are good enough up 
too bombarding energies of about 10 AGeV. Due to the projectile-target transparency, they 
become less and less justified with increasing bombarding energy. To take into account 
this effect, generalized multi-fluid models have been constructed in Refs. 
important ingredients of such models, the inter-fiuid coupling terms, are rather uncertain 
and usually parameterized phenomenologically. 

It is clear that hydrodynamical approach can not be directly applied to late stages of a 
heavy-ion reaction when binary collisions of particles become too rare to maintain the local 
thermodynamic equilibrium. A standard way 10| to circumvent this difficulty is to introduce 
a so called "freeze-out" criterion to stop the hydrodynamical description. Often it is pos- 
tulated that collisionless expansion of particles starts at some isothermal hypersurface 41|. 
Unfortunately, this approximation is rather crude [42] and even in contradiction with exp er- 



imental data (see e.g. 



43|). A more consistent procedure has been suggested 
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46| 



within a hybrid "hydro-cascade" model. In this scheme hydrodynamics is used for generat- 
ing coordinates and momenta of hadrons at some intermediate stage of the reaction. These 
characteristics are used for transport simulations of later stages. 

Below we formulate a version of the ideal (3+l)-dimensional hydrodynamics suitable 
for the domain of AGS, NICA, FAIR and SPS energies. This model belongs to the first 
class and uses a new EoS {4?! with the deconfinement and liquid-gas phase transitions. 
In our calculations we perform a detailed analysis of the dynamics of nuclear collisions at 
various energies giving a particular attention to macroscopic characteristics, i.e. collective 
flows, life-time and volume of quark-gluon and mixed phases. It is shown that maximal 
energy and baryon densities predicted by our fluid-dynamical simulations agree well with 
ID shock wave calculations. A special analysis is carried out to evaluate the importance 
of transparency effects. This is made by comparing our results with predictions of the 
three-fluid model of Ref. j4oj]. To calculate momentum spectra as well as directed and 
elliptic flows of secondary particles, we apply an approximation of isochronous freeze-out. 
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In these calculations we investigate the sensitivity of observables to finite size of hadrons by 
introducing the excluded volume corrections. 

The paper is organized as follows: our hydrodynamic model is formulated in Sec. II. The 
results of numerical calculations for central Au+Au collisions at FAIR and SPS energies are 
given in Sec. Ill . Here we also compare the predictions of the one- and three-fiuid models. 
In Sec. IV we discuss a space-time picture of non-central Au+Au collisions. Particle spectra 
and parameters of transverse collective fiows are considered in Sec. V and VI. Our conclusions 
are presented in Sec. VII. 



II. FORMULATION OF THE MODEL 



A. Equations of ideal fluid dynamics 



Below we study the evolution of highly excited, and possibly deconfined, strongly- 
interacting matter produced in ultrarelativistic heavy-ion collisions. It is assumed that 
this evolution can be described by the equations of ideal relativistic hydrodynamics 48|. 
These equations represent local conservation laws of the 4-momentum and baryon charge 

5,T^'^ = 0, (1) 
d^N'' = 0. (2) 

In the limit of small dissipation the baryon 4-current A^'^ and the energy-momentum ten- 
sor T^^ can be expressed as {h = c = 1) 

= nu^, (3) 

T^" = {e + P)u^'u'' -Pg"", (4) 

where e, n and P are the rest-frame energy density, the net baryon density and pressure of 
the fiuid, = 7(1, v)^ is its collective 4-velocity, 7 = (1 — v^)"^/^, gf^" = diag (+, — , — , — ) 
is the metric tensor. Here and below we denote by v the fiuid 3-velocity. 
One can rewrite Eqs. (II])-® in Cartesian coordinates (t, r) as follows 

dtE +V{vE) = -V{vP), (5) 
dtM +V{vM) = -VP, (6) 
dtN+V{vN) = 0, (7) 



where E = T^^, M'^ = T^^ and = are the energy density, the 3- momentum density 
and the baryon density in the lab frame. 

The relations between hydrodynamical variables in the lab and local rest frames can be 
written as 

E = ^^{e + v''P), (8) 
M = -i'^{e + P)v, (9) 
N = -in. (10) 

To solve Eqs. fl5|)- f|T0|) it is necessary to specify the EoS P = P{n,e) of the fluid and initial 
conditions. 



B. Equation of state 

In our calculations we use the EoS of strongly interacting matter with the first order 



deconfinement phase transition (EoS-PT) 47|]. The hadronic phase (HP) is described as 
the hadron resonance gas with inclusion of known hadrons with masses up to 2 GeV. Finite 
size effects are taken into account by introducing the excluded volume corrections. The 
same excluded volume parameter f e = 1 fm^ is used for all hadronic species. A Skyrme-like 
mean field for baryons 

U{n) = -an + I3n^/'^ (11) 
is added to guarantee that the hadronic matter has correct saturation point and a liquid-gas 



phase transition 72] ■ The deconfined quark-gluon phase (QGP) is described by the bag 
model EoS with lowest-order perturbative corrections. The phase transition boundaries and 
characteristics of the mixed phase (MP) are found by using the Gibbs conditions with the 



strangeness neutrality constraint (for details, see Ref. 47|]). The domains of MP states in 
different thermodynamic variables are shown in Figs. [THU [6H71 [9]-[T0l To probe sensitivity to 
the EoS, we have also performed calculations with the EoS of ideal hadron gas (EoS-HG). 
In this case we disregard the excluded volume effects assuming that fg = 0. To provide 
stability of initial nuclei we also introduce the Skyrme-like mean field. The parameters of 
EoS-PT and EoS-HG are given in jlj. 



The graphic representation of P{n,e) for both EoSs is given in Figs. [TH2] 73|. One can 



see that values of pressure predicted by these two EoSs significantly differ at large enough e 




FIG. 1: (Color online) Pressure as a function of energy density at different baryon densities, scaled 
by saturation nuclear density uq = 0.15 fm^^. Solid and dashed lines correspond to EoS-PT and 
EoS-HG, respectively. The MP region of the deconfinement phase transition is shown by shading. 

and n. An important information about the dynamics of matter in the heavy-ion colhsion 
can be drawn from Fig. [2] which shows adiabatic trajectories, i.e. hnes of constant entropy 



per baryon, a = s/n, where s is the entropy density. As well known 48|], this quantity is 
conserved along the trajectories of fluid elements in the ideal hydrodynamics (in absence of 
shock waves). Our calculations show that larger values of a in central heavy-ion collisions 
are achieved at larger bombarding energies -Eiab- The slopes of adiabates in Fig. [2] give the 



sound velocities squared 



48| 



From from Fig. [2]one can clearly see that Cs reach minimal values (" softest points" ) in the MP 
region. This should essentially influence the expansion dynamics. Indeed, an instantaneous 
acceleration of a fluid element is proportional to the gradient of pressure VP = c^Ve . At 
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FIG. 2: (Color online) Adiabatic trajectories in the e — P plane. Thick and thin lines correspond 
to EoS-PT and EoS-HG, respectively. Shading shows the MP region of the deconfinement phase 
transition. 

given Ve one may expect smaller accelerations of matter for "softer" EoSs, which have 
smaller values of sound velocity. From Fig. [2] one can see that although generally the 
EoS with the phase transition is softer than the EoS of pure hadron gas, in some regions of 
thermodynamic parameters the EoS-PT has sound velocities of the same order or even larger 
than those in the EoS-HG. Such a situation may occur in the MP region for intermediate 
values of specific entropy. This anomaly is expected in the NICA-FAIR energy domain, 
where our calculations predict a non-monotonous behavior of collective flows and particle 
spectra as functions of bombarding energy (see below). 

C. Initial conditions 

We start our hydrodynamical simulation from the stage when two cold nuclei approach 
each other. Unless stated otherwise, we consider collisions of gold nuclei {Z = 79,A = 197) 
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which have the Woods-Saxon density distribution in their rest frame, 

^(^■,^=0) = , , (13) 

1 + exp (^^^ 

with the parameters hq = 0.15 fm~^, a = 0.5 fm, Rq = 6.7 fm. We assume that the rest 
frame initial energy density of each nucleus, e{r,t = 0) is given by the r.h.s. of Eq. f[T^ with 
the replacement of uq by Eq = fiQUQ where /io = 0.923 GeV is the baryon chemical potential 
of equilibrium nuclear matter. Due to the stabilizing effect of the Skyrme-like mean field, 
the initial nuclei may stay in equilibrium with vacuum at P = . They propagate without 
distortion until their density distributions essentially overlap. 



The cm. collective velocity of nuclei at t = is taken as Vq = A/-Eiab/(2m7v + -Eiab), 
where ttln = 0.939 GeV is the nucleon mass and -Eiab is the projectile bombarding energy 
(per nucleon) in the lab frame. Below the beam axis is denoted by z and the x-axis is 
chosen along the impact parameter vector b. In this case the reaction and transverse planes 
correspond to y = and z = 0, respectively. The initial longitudinal distance between the 
target and projectile centers in the cm. frame is chosen as 

zt(0)-^p(0) = 2(i?o + 6.9a)/7o, (14) 
where 70 = (1 — fo)~^^^ is the Lorentz factor of colliding nuclei j?^. 

D. Numerical scheme 

The numerical solution of fluid-dynamical equations (l5])-(I7j) was obtained by using the 



flux-corrected transport algorithm j49|, |50|. All calculations are performed in the cm. frame. 
Unless otherwise stated, we use a cubic Eulerian grid with the cell size Aa; = 0.1 fm along 
each direction. Typically we take the numbers of grid points equal to 500 x 300 x 400 in 
X, y, z directions and choose the time step At = 0.01 fm/c. In our simulations we use linear 
interpolations of tables P(n, e) prepared with fixed steps in n and e. We have checked that 
our numerical scheme conserves the total baryon number and energy with relative accuracy 
better than 1%. The program code explicitly uses the symmetry of the collision process with 
respect to the reaction plane y = 0. This reduces the computational memory by a factor 
of two. A typical run for the 10 AGeV Au+Au collision at 6 = 7 fm requires about 6 Gb 
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FIG. 3: (Color online) The rest-frame baryon density n as a function of time in a central (b = 0) 
Au+Au collision at -Eiab = 10 AGeV. Dotted lines show density in a central cell, while solid ones 
represent density values averaged over the central box. Different colors correspond to calculations 
with different Ax. All results are obtained for EoS-PT. 

of memory and approximately 6 hours of CPU time on a machine with the 2.1 GHz AMD 
Opteron processor. 

Below we often analyze time evolution of fluid-dynamical quantities averaged over a 
"central box" around the symmetry point x = y = z = 0. As a rule, we use the box with 
dimensions 75 1 

\x\, \y\, 7o|z| < Ifm. (15) 

For comparison Fig. [3] shows the time dependence of the baryon density in the central grid 
cell and in the central box for the case of a central Au+Au collision at i^iab = 10 AGeV. 
One can see that using rough grids leads to oscillations of the energy density which are of 
numerical origin. They are especially visible at early stages when a shock-like compression 
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of matter takes place (see Sec. IIIB). Averaging over neighboring cells partly remove these 
oscillations. 



E. Particle spectra and parameters of collective flows 

To calculate hadronic momentum distributions we apply the approximation of instan- 
taneous freeze-out: it is assumed that a sudden transition from the local equilibrium to 



76|. 



coUaonte. propagation of particle, take, place at .on,e .pac^tin,e hype^urface 
Within this approximation one can use the standard Cooper-Frye formula [41] for the in- 
variant momentum distribution of the hadronic species i 

^3^ = X3^ = 77^ / ^^M^M^xp ^ ±U , (16) 



d^p dyd^pT (27r)3 7 \ T 

where p'^ = {E, p)^ is the 4-momentum of the particle, y and Pt are, respectively, its lon- 
gitudinal rapidity and transverse momentum, T denotes the local temperature, Qi is the 
statistical weight of the hadron species i. 

Plus and minus in Eq. ( JT6l) correspond, respectively, to fermions and bosons. Using con- 
ditions of chemical equilibrium one can express the particle's chemical potential Hi through 



the baryon (/x) and strange (/x^) chemical potentials as follows |47l. |51| 



Hi = BiH + SiPs , (17) 

where Bi and Si are, respectively, the baryon and strangeness number of species i . 

In calculating spectra of baryons {Bi = 1) we take into account the mean field effects. 
In this case we use Eqs. (IT6l)-([T7l) with the replacement of the baryon chemical potential h 
by its "kinetic" part fiK = A* ~ U{n) where U{n) is the mean field potential introduced 
in Eq. f fTTl) . As indicated above, our calculations with the deconfinement phase transition 

n 

use the EoS with inclusion of finite volume corrections. In accordance with Ref. [47|, we 
introduce these corrections by the replacement /ij — )■ /ij — ViPx where is the kinetic 
part [47] of the total pressure and Vi is the excluded volume of i-th hadron. Unless stated 
otherwise, particle spectra and collective flows in the case of EoS-PT are calculated by using 
the parameter Vi = 1 fm'^ for all hadronic species. 

For our qualitative analysis we choose the simplest condition of isochronous freeze-out 
(t = tfr = const). Then da^ = d^x6^fi and Eq. (fT6|) gives the hadronic momentum dis- 
tribution at a fixed time t in the cm. frame. In addition to contributions of "thermal" 
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nucleons and pions, which are calculated directly by using Eq. (fT6|) . we also take into ac- 
count resonance decays, e.g. A — )■ Ntt, p — t- 27r. Below we assume zero widths of resonances 



and apply t 
R-^i+X 



le following formula 16[ for the spectrum of hadrons i from two-body decays 
77| 



E- 



d^PR 



d^N, 



R ^ I PPr 



Ef 



(18) 



d^p Anqo J d^pR \ tur 

Here rriR and pr are the mass and 4-momentum of the resonance R {R = A, p . . .) , Eq and go 

are the energy and 3-momentum of the hadron i in the rest frame of R, hi is the branching 

ratio of the decay channel R-^i+X. Equation f|T6|) is used to calculate spectrum of 

resonances (in the Boltzmann approximation). All known hadronic resonances with masses 

up to 2 GeV are taken into account. The statistical weights, branching ratios and masses of 

these resonances are taken from 52 1. 



Transverse collective flows of matter created in heavy-ion collisions, are rather sensitive 
to its EoS. Much attention in recent years is given to the parameter of elliptic flow V2 
The 2+1 fluid dynamical calculations at SPS and RHIC energies 
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53|. 



54| show a significant 



sensitivity of this quantity to the deconfinement phase transition. The elliptic flow of i-th 
hadrons {i = tt,N . . .) is usually determined as 



vf{y) 



d^PT cos (24)) Ed^Ni/d^p 



d^PrEd^Ni/d^p 



(19) 



where = aiccos {p^/pt) is the azimuthal angle of px with respect to the reaction plane. 
Here d^Ni/d^p is the momentum distribution of i-th hadrons with inclusion of resonance 
decays. Unfortunately, the elliptic flow is sensitive 45I, 55, 56| not only to the EoS, but also to 
dissipative and freeze-out effects. Below we also calculate the directed flow parameter Vi\y). 
It is defined by the r.h.s. of Eq. (IT^ with the replacement cos (20) — )■ cos 0. It is clear that 
f 1, f2 — )■ for purely central collisions (6 — > 0) . 

To discuss qualitatively possible differences between the EoS-PT and EoS-HG at different 
bombarding energies, below we calculate a so-called momentum anisotropy parameter 
which characterizes the collective flow asymmetry in the transverse plane. Following Ref. 
we define this quantity as 

''dxdy (T^^ -T^^) 

(20) 



dxdy {T^^^T^y) 
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where T^^, T^^ are the components of momentum density in the transverse plane z = 78 |. 



In Ref. the approximate relation ep ~ 2^2'^^ (0) has been obtained from (2+1)- 

hydrodynamical simulations of Pb+Pb collisions at SPS and RHIC energies. 



III. EVOLUTION OF FLUID DYNAMIC QUANTITIES IN CENTRAL COLLI- 
SIONS 



A. Comparison of results for EOS-PT and EOS-HG 




FIG. 4: (Color online) Time evolution of the energy density in the central box of a central Au+Au 
collision for various bombarding energies. Solid and dashed lines correspond to EoS-PT and 
EoS-HG, respectively. 

Below we discuss thermodynamic characteristics of matter created in Au+Au collisions 
at various bombarding energies. Our goal is to find differences between the results obtained 
with EoS-HG and EoS-PT. Let us first consider purely central collisions with 6 = 0. Fig- 
ures HHH show the time evolution of the rest frame energy and baryon densities averaged over 
the central box defined in Sec. Ill D[ One can see that after a rapid initial growth the baryon 
and energy densities remain practically constant during certain time intervals (depending on 
the bombarding energy -Eiab) and only later they start to decrease. The plateau is explained 
by running shock waves created at an early stage of the reaction. Within the considered 
one-fluid approach the entropy of matter is generated via the shock-like compression mech- 
anism (see the next section) which converts the cm. kinetic energy of projectile and target 
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nucleons into internal energy of stopped matter. The plateaus in Figs. |1H5] begin at the 
time moments when shock waves reach dense central parts of colliding nuclei. Initially these 
shocks appear at the central point r = in the cm. frame and then propagate outwards 
through the colliding nuclei. Right-end points of the plateaus correspond to the moments 
when rarefaction waves, propagating from rear sides of nuclei, reach the central point. One 
can see that at i?iab ^ 5 AGeV, when matter enters MP states, the calculation with the 
EoS-PT predict larger values of e and n as compared with the EoS-HG. 



1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — |- 
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2 4 6 
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FIG. 5: (Color online) Same as Fig. HJ but for net baryon density. 



B. Simplified picture of matter evolution in a central box 



To estimate parameters of states with maximal compression, let us consider a collision 
of two slabs of cold nuclear matter (n = uq^T = {]) moving initially with velocities ±fo 
along the z-axis. After the first contact two shock fronts start to propagate with constant 
velocities in the positive and negative z-directions. From the continuity of T^^ ,T^^ , in 
the shock front rest frame one gets 48|] the relation (so-called " Taub adiabate" ) connecting 
the baryon (n) and energy [e) densities of fluid behind the shock front: 



' 



(21) 



where = A^o'^o is the initial rest-frame energy density of each slab [79 1. 
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FIG. 6: (Color online) Taub adiabates and parameters of shock waves for collisions of two slabs 
of normal nuclear matter for -Eiab = 10 and 40 AGeV. Thin straight lines represent stopping 
conditions ()22p for these values of bombarding energy. Thick solid and dashed lines are the Taub 
adiabates (j2ip for EoS-PT and EoS-HG, respectively. The full dot shows the initial state n = 
nQ,T = 0. Open dots represent parameters of compressed zone behind the shock fronts. Shading 
shows the MP region of deconfinement phase transition. 

In addition one can write down the equation expressing conservation of energy per baryon 
at the shock front. In the cm. frame this equation looks as the stopping condition 

e/n = 7oeo/^o, (22) 

where 70 = (1 + -E'lab/SmAr)^/^ is the initial Lorentz factor. Graphic solutions of 
Eqs. fl2T]) -f l22|) for two EoSs used in this paper are shown in Fig. |6l One can see that 
at -Eiab ^ 10 AGeV calculations with the EoS-PT predict noticeably higher maximal values 
of e and n as compared with the EoS-HG. This in turn leads to higher collective flows 
of matter in hydro simulations with the phase transition. Comparison with parameters of 
maximal compression in central Au+Au collisions (see Figs. HHS]) shows a good agreement 
between the Taub adiabate predictions and fluid dynamical calculations. Our shock wave 
calculations show that phase transition effects and, in particular, excitation of MP states 
may be expected only at sufficiently high bombarding energies, i?iab ^ 3 AGeV [80,] • 
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12 3 



e (GeV/fm^) 

FIG. 7: (Color online) Pressure as a function of energy density in a central box for central 10 AGeV 
Au+Au collisions. Dashed and solid lines give the results for EoS-HG and EoS-PT, respectively. 
Numbers in squares and circles denote the values of cm. time (shown are only states at i > 6 fm/c). 
Shading represents the MP region. 

Decreasing parts of curves in Figs. HHS] correspond to the expansion stage of system evo- 
lution in the central box. It can be shown that this expansion proceeds quasi-adiabatically, 
i.e. with nearly constant entropy per baryon a (small deviations are caused by numerical 
viscosity and by averaging over cells in the central box). Figure [7] shows time evolution of 
thermodynamic parameters in the central box for the case of central Au+Au collisions at 
E\ah = 10 AGeV. Only late stages of the reaction in the e — P plane are considered. Com- 
parison with trajectories in Fig. [2] shows that in the considered case the fluid expansion in 
the central box proceeds along the adiabate with a ~ 10 (for both EoSs). From Figs. HHSllT] 
one can see that in the case of the EoS-PT the system expands with some delay as compared 
with the EoS-HG, i.e. the same values of e and n are achieved at later times. However, 
due to differences in pressure, final hadronic states predicted for the EoS-PT and EoS-HG 
have, generally speaking, different collective velocities. 

C. Comparison of one and three fluid models 

It should be stressed that the results presented above were obtained under the assumption 
of full mutual stopping of projectile and target matter. Presumably, this approximation 
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15 AGeV 
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FIG. 8: (Color online) Time evolution of energy density in the central box of central Au+Au 
collisions with different i?iab- Solid lines shows our results in the case of EoS-HG and sharp 
density distributions of initial nuclei. Dashed lines represent three-fluid calculations [59 1. 



should be less justified with raising bombarding energies (see a more detailed discussion 
in js^]). To estimate the importance of projectile-target transparency effects, we have 



compared our calculations with the results of the three-fluid model of Ref . [40| . The standard 
version of this model applies a purely hadronic EoS. In addition, sharp density profiles of 
initial nuclei are assumed. To make the comparison possible, we performed our calculations 
with the EoS-HG choosing a small diffuseness parameter a = 0.01 fm (see Eq. f|T3|) ). The 
results of this comparison are presented in Fig. [HI One can see that in the case of central 
Au+Au collisions the differences between the one- and three-fluid calculations are rather 
small at -Eiab ^ 30 AGeV. At larger bombarding energies maximal values of e predicted by 
one-fluid calculations are about 10% larger than in three-fluid model. In addition, the latter 
predicts a noticeably slower increase of energy density at the initial stage of the collision. 
From these results we conclude that transparency effects are not so strong at least in central 
Au+Au collisions at FAIR energies and below. 



D. Space time picture of central Au+Au collision 

Figures [OHTOl represent dynamical trajectories of matter (the central box) produced in 
central Au+Au collisions at different bombarding energies. Shown are results for n — e 
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FIG. 9: (Color online) Time evolution of matter in central Au+Au collision at different -Biab- 
Shown are values of energy and baryon densities averaged over the central box. Dashed and solid 
lines correspond to EoS-HG and EoS-PT, respectively. Numbers in triangles and circles give the 
cm. time in fm/c. Shading shows MP region of deconfinement phase transition. 

and yU — T planes, respectively. As compared with the EoS-HG, the calculations with the 
phase transition predict longer life-times of states with maximal compression and additional 
time delays of expansion stages. Especially large differences between two EoSs are visible at 
intermediate times in the fi — T plane. The dotted line in Fig. [TOl shows the MP region in this 



plane 8l|. The calculations with EoS-PT predict a zigzag-like behavior of the trajectories 
in the fi — T plane with slightly raising temperature in the MP as a function of time. On the 
other hand, the trajectories do not bend at the MP boundaries in the n — e plane. According 
to Figs. lOHTOl the largest life times of the MP in central Au-|-Au collisions, about 4 fm/c, 
are reached at E'lab ~ 10 AGeV. At lower energies the MP region is traversed only partly, 
while at higher Ei^y, the produced matter expands too fast. 

Figures [TTHT2] shows the temperature and pressure profiles as functions of the longitudinal 
coordinate along the beam axis {x = y = 0) for central Au+Au collision at i?iab = 10 AGeV. 
It is interesting that in the case of EoS-PT the temperature at small \z\ changes non- 
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FIG. 10: (Color online) Time evolution of matter (central cell) in the l-i — T plane. Left and right 
panels correspond to central Au+Au collisions at i^'iab = 10 and 20 AGeV, respectively. Dashed 
(solid) lines show the results for EoS-HG (EoS-PT). Numbers in squares and circles give time 
values in fm/c. 
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FIG. 11: (Color online) The temperature profiles along the beam axis at different time moments 
in a central 10 AGeV Au+Au collision. The left (right) panel corresponds to EoS-PT (EOS-HG). 
The results for z > are shown only. 

monotonically in the interval t = 8 — 12 fm/c. This follows from a specific behavior of the 
critical temperature as a function of fi. Indeed, as one can see from Fig. [10], the temperature 
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FIG. 12: (Color online) Same as Fig. [TT] but for the pressure profiles along the beam axis. 




FIG. 13: (Color online) Contours of equal energy densities (in GeV/fm^) in the reaction {y = 0) 
plane at t = 8 fm/c for central 10 AGeV Au+Au collision. The left (right) panel corresponds to 
EoS-PT (EoS-HG). Numbers near the vertical and horizontal lines show the x and z coordinates 
in fm. Thick line in the left panel shows the boundary of MP. 

of MP in a central slice increases noticeably during the expansion. Such a behavior is 
explained by the release of the latent heat in hadronizing process. According to Fig. [121 
the deconfinement phase transition leads to slower drop of pressure gradients in expanding 
matter along the beam direction. These gradients are especially large at the MP boundary. 
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This results in strong acceleration of fluid cells during the hadronization processes. As 
will be shown below, such acceleration leads to additional broadening of baryon rapidity 
distribution as compared to a purely hadronic scenario (see Sec. V). 

Comparison of spatial distributions of energy density for two EoSs is made in Fig. [121 It 
reprents the contour plots of e in the reaction plane of a 10 AGeV central Au+Au collision 
at t = 8 fm/c. The calculation with the EoS-PT shows that at this time the central zone 
is still in the MP. In this case maximal values of energy density are noticeably higher as 
compared with the EoS-HG. 

IV. MODELING NON-CENTRAL COLLISIONS 

A. Energy density and collective velocities in Au+Au collisions at E'lab = 10 AGeV 

Let us consider now non-central collisions. As we shall see below, especially interesting 
is the region of bombarding energies ~ 10 AGeV. At such energies we expect an enhanced 
sensitivity of collective flow observables and particle spectra to the phase transition. In 




FIG. 14: (Color online) Contours of equal energy densities in the reaction [y = 0) plane at 
t = 8 fm/c for 10 AGeV Au+Au collision with 6 = 7 fm. The left (right) panel corresponds to 
EoS-PT (EoS-HG). Numbers near the vertical and horizontal lines show the x and z coordinates 
in fm. Thick line in the left panel shows the boundary of MP. 
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FIG. 15: (Color online) Same as Fig. [T4l but for the transverse plane z = 0. 
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FIG. 16: (Color online) Collective velocities of matter in the reaction plane for midcentral 10 AGeV 
Au+Au collision (t = 8 fm/c, 6 = 7 fm). Left and right panels correspond to EoS-PT and EoS-HG 
respectively. Lengths and directions of arrows represent 3-velocity vectors. Numbers show the z 
and X coordinates in fm. Thick line in the left panel shows the boundary of MP. Vertical color 
strips characterize modules of 3-velocity in units of c. 

this section we present results for 10 AGeV Au+Au collisions with the impact parame- 
ter 6 = 7 fm. 

Figures [HHTT] show contour plots of the energy density and 3- velocity fields in the reac- 
tion {z — x) and transverse {y — x) planes at t = 8 fm/c. One can see that in calculation 
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with the EoS-PT a significant part of the system at this stage is still in the MP. 




FIG. 17: (Color online) Collective velocities of matter in the transverse plane for 10 AGeV Au+Au 
collision (t = 8 fm/c, 6 = 7 fm). Upper and lower panels correspond to EoS-PT and EoS-HG 
respectively. Numbers show the x and y coordinates in fm. Thick line in the upper panel shows 
the boundary of MP. Vertical color strips characterize modules of 3-velocity in units of c. 

Figures [TBHT71 represent the collective velocity fields for the same reaction. Blue colored 
central regions in these figures correspond to spatial domains with relatively low absolute 
values of cm. velocity v < 0.2 c. At the considered time moment, such a region occupies 
much larger fraction of space in the calculation with the deconfinement phase transition. Ac- 
cording to Figs. [T6HT71 the MP matter is practically at rest in the cm. frame. On the other 
hand, pronounced collective flows of hadronic matter are formed in transverse directions. 
Especially large transverse velocity components appear in the region of geometrical overlap 
of initial nuclei, at outer edges of the HP. Qualitative structure of the velocity flelds may 
be understood if one takes into account that fluid elements are accelerated mainly along the 
directions of largest gradients of energy density (see Sec. IIIB|) . In particular, this leads to 



the effect of "antiflow" i.e. to appearance of fluid elements with opposite signs of Vx and 
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FIG. 18: (Color online) The pressure profiles in 10 AGeV Au+Au collision (6 = 7 fm). Left and 
right panels correspond, respectively, to the x and y axis in the transverse plane z = 0. Red 
and blue curves give the results for t = 8 and 9 fm/c. Solid and dashed lines are calculated with 
EoS-PT and EoS-HG, respectively. 

According to Fig. [151 at i^iab — 10 AGeV the EoS-PT leads to an essentially different 
spatial symmetry of the energy density distribution in the transverse plane as compared 
with the purely hadronic EoS. On the other hand, the calculation with the EoS-PT predicts 
much larger pressure gradients and accelerations of matter in the x and y directions. In 
Fig. [18] we compare the pressure profiles for both EoSs at different times. One can see, that 
in the case of the EoS-PT, especially strong pressure gradients at t ~ 8 fm/c are formed 
in the x direction, at the boundary between the hadronic and mixed phases. This in turn 
leads to an enhanced momentum anisotropy in the case of EoS-PT as compared with the 



EoS-HG (see next section) [82| 



B. The momentum anisotropy at different bombarding energies 

Here we present the results concerning the momentum anisotropy in (semi)peripheral 
Au+Au collisions. This quantity has been defined in Eq. (j20l) . Figure [TOl represents the 
time evolution of Cp at different i^iab • Thick lines are calculated using the generalized 
definition of with y^, = 0.5 (see footnote [7]). By thin lines (for E'lab = 10 AGeV) we also 
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present calculations in the limit — )■ 0. One can see that in the case with phase transition 
the simplified formulae f l20|) leads to unphysical maximum of ep{t) at intermediate times. 
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FIG. 19: (Color online) Time dependence of the momentum anisotropy in Au+Au collisions at dif- 
ferent bombarding energies (6 = 7 fm). The solid (dashed) lines correspond to EoS-PT (EoS-HG). 
Dots marks time moments when the energy density in central box becomes lower certain values e^. 
Thin lines in the left panel show the results of calculations for -Eiab = 10 AGeV with simplified 
definition of ep (see text). 

According to Fig. [191 the momentum anisotropy is rather sensitive to EoS in the bombarding 
energy range -Eiab ^ 20 AGeV. The asymptotic values of ep at -Eiab ^ 10 AGeV are larger in 
calculations with the deconfinement phase transition. We think that such "counterintuitive" 



result 83| may be explained by the difference in pressure gradients in the transverse plane. 
As discussed in the end of preceding section, they are much larger at intermediate times if 
the MP is created. 

As mentioned in Sec. Ill El the freeze-out effects may significantly change the observed 
values of transverse fiows in heavy-ion collisions. To estimate a possible spread of results we 
determine the time moment tf^ when the energy density in a central box becomes smaller 
a certain freeze-out value 5^. Applying such freeze-out condition, in Fig. [19] we mark by 
points the values of corresponding to different values of e^. Figure [20] shows excitation 
functions of the momentum anisotropy for two EoSs. The lines connect the values of ep 
taken at t = where corresponds to = 0.4 GeV/im^. At each i^iab the ends of arrows 
in Fig. [20l show the values of ep for = 0.2 and 0.6 GeV/fm^. In the case of EoS-PT our 
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FIG. 20: (Color online) Excitation function of momentum anisotropy in Au+Au collisions with 
6 = 7 fm. The solid (dashed) lines correspond to EoS-PT (EoS-HG). Arrows show possible shifts 
of ep-values for different choices of freeze-out energy density between 0.2 and 0.6 GeV/fm'^. 

model predicts a non-monotonic dependence ep(£'iab) with maximum at i?iab — 10 AGeV. 

As indicated above, the momentum anisotropy is approximately proportional to the pion 
elliptic flow at midrapidity. It is interesting that existing experimental data on the proton 
and pion elliptic flows in Au+Au and Pb+Pb collisions (see their compilation in Ref. jssl) do 
not exclude presence of a local maximum of V2 at high AGS energies. However, the data at 
AGS and SPS energies were obtained by using different detectors, centrality cuts and meth- 
ods of f 2 determination. It would be highly desirable to perform more detailed measurements 
of the elliptic flow excitation function in the bombarding energy range 5-50 AGeV. Hope- 
fully, low energy RHIC as well as future NICA and FAIR experiments will help in obtaining 
such data. 

V. HADRONIC SPECTRA 

In this section we present the model results for proton and pion momentum distribu- 
tions in central Au+Au collisions at i^iab — 10 A GeV. These distributions are calculated 
by using the Cooper-Frye formula (fT6|) with the isochronous freeze-out hypersurface. The 
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FIG. 21: (Color online) Transverse kinetic energy distributions of protons (left panel) and negative 
pions (right panel) in 10.7 AGeV Au+Au collision (6 = 2 fm). Thick solid and dashed lines 
are calculated with EOS-PT and EOS-HG, respectively. Thin lines show the distributions of 
thermal hadrons, i.e. without contributions from decays of resonances. Experimental data are 



taken from 



6l|. 



results for the EoS-PT include the excluded volume corrections (see Sect. Ill Ep . The proton 
and 7r~ distributions are obtained from nucleon and pion spectra by introducing the scaling 
factors 1/2 and 1/3, respectively. The parameter of freeze-out time tfr is chosen to achieve 



the best fit of experimental data [8J|. Our calculations with the EoS-PT show that the best 
choice corresponds roughly to the end of hadronization stage. In this case the thermody- 
namic parameters of central cells become close to the boundary between the MP and HP 
(see Figs.fflinD. 

According to our analysis, significant fractions of proton and pion yields in the considered 
reaction come from decays of resonances. At -Eiab ~ 10 AGeV, relative contributions of free 
("thermal") pions and protons at midrapidity are about 25 and 60 percents of total yields 
of pions and protons, respectively. Especially important are contributions of isobar decays. 
This is explained by large values of baryon chemical potential (yU ~ 0.6 — 0.7 GeV at freeze- 
out) expected in nuclear collisions at AGS and FAIR energies. Having in mind that measured 
yields include feeding from weak decays of hyperons and mesons, we take into account 
the corresponding contributions in our calculations. 
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FIG. 22: (Color online) Rapidity distributions of protons in 10.7 AGeV Au+Au collision (6 = 2 fm). 



Solid and dashed 
experimental data 



ines are calculated with EOS-PT and EOS-HG, respectively. Full dots are 
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631 ] ■ circles are obtained by reflection with respect to midrapidity. 



According to Fig. EH the proton and pion pt distributions in central Au+Au collisions at 
Eiah ~ 10 AGeV are not so sensitive to EoS. In this figure we explicitly show the contribu- 
tions of " directly" produced (thermal) proton and pions. Indeed, one can see that resonance 
decays give large contributions at all px, especially for pions. 

The proton rapidity distribution for the same reaction is presented in Fig. |22l As com- 
pared to the EoS-HG, the calculation with the phase transition predicts a noticeably broader 
rapidity distribution. In this case the agreement with experimental data outside the central 
rapidity region is better than in the purely hadronic scenario. The physical meaning for this 
broadening is related to larger pressure gradients in the longitudinal direction, as demon- 
strated in Fig. [121 Our analysis shows that the shape of the proton rapidity distribution at 
small cm. rapidities is rather sensitive to the excluded volume parameter. It is interesting 
that single- and three-fluid calculations of central 11 AGeV Au+Au collisions also pre- 
dict appearance of a deep in the proton rapidity distribution at midrapidity for a EoS with 
the deconfinement phase transition. Irregular behavior of this distribution as a function of 
bombarding energy has been recently suggested {gS] as a signature of the QGP formation. 
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FIG. 23: (Color online) Same as Fig. [20] but for rapidity distributions of vr" mesons. Full sym- 
bols are experimental data for charged pions in central Au+Au collisions at Slab = 8 [66] and 
9.9 AGeV. Open symbols are obtained by reflection. 

The TT" rapidity distribution for the 10.7 AGeV central Au+Au coUision is shown in 
Fig. [231 In the case of the EoS-PT, a reasonable agreement with experimental data can be 



achieved only if we assume a freeze-out time smaller than for protons [85[. Again, calculations 
with the phase transition predict a broader distribution than those for the EoS-HG. However, 
differences in shapes of pion rapidity distributions are not so strong as for protons. It is 
worth to note that our results on hadronic spectra and collective flows (see next section) 
are rather preliminary. First, they are obtained by using the simplest option of isochronous 
freeze-out. Second, the resonance decays are considered rather schematically i.e. the width 
of resonances was neglected and multi-particles decays were included only approximately. 



VI. COLLECTIVE FLOWS 

The momentum anisotropy considered in Sec. ilVBl is not measurable directly and can be 
used only for qualitative analysis of collective flows in nuclear collisions. On the other hand, 
theoretical and experimental studies of the directed (fi)and elliptic (^2) flows became very 
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FIG. 24: (Color online) Directed flows of protons (left panel) and charged pions (right panel) as 
functions of cm. rapidity in 10 AGeV Au+Au collisions. Solid (dashed) lines are calculated with 
EoS-PT (EoS-HG). Thick and thin lines correspond to 6 = 4 and 7 fm, respectively. 



popular in recent years (see e.g. |67|). In particular, this was initiated by the prediction [33 1 
of a "softest point collapse" of the nucleon directed flow at E'lab ~ 8 AGeV. In this section 
we present our results for Au+Au and Pb+Pb collisions at -Eiab = 10 and 40 AGeV. These 
results are obtained by using formulas of Sec. IIIEi 

Figure EH shows the model predictions for the proton and pion directed flows in 10 AGeV 
Au+Au collisions. As before, we compare the results obtained for the EoS-PT and EoS-HG. 
To demonstrate sensitivity to the centrality, we present results for two values of the impact 
parameter, 6 = 4 and 7 fm. In both cases it is assumed that freeze-out occurs when the 
energy density in the central cell becomes lower than Ef^ = 0.3 GeV/fm^. The model predicts 
a strong reduction of |fi | at small \y\ in transition from the EoS-HG to EoS-PT. One can 
see a noticeable difference of slopes dvi/dy between two EoSs in a central rapidity region. In 
addition, negative slopes (antifiow) of the pion and proton vi are predicted in the calculation 
with the EOS-PT. However, this effect is stronger for pions. Our analysis shows that the 
antifiow in this case is formed mainly due to the contribution of thermal pions. Note, that 
slopes dv^^/dy at \y\ < 0.5 only slightly depend on b. 

Unfortunately, the rapidity dependence of directed flows is not measured at AGS bom- 
barding energies. On the other hand, there are measurements [gsI of a similar quantity, 
the so-called "sideflow" <px>, which is defined as a mean transverse momentum (per 
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FIG. 25: (Color online) Mean directed transverse momentum of protons as functions of cm. 
rapidity in 10 AGeV Au+Au collisions (6 = 4 fm). Solid and dashed lines are calculated with 
EoS-PT and EoS-HG, respectively. Experimental data [68] for protons (full dots) and nucleons 
(full triangles) correspond to 10.1 AGeV midcentral Au+Au collisions. Circles are obtained by 
reflection. 

particle) projected on the reaction plane. Below we calculate this quantity by using the 
r.h.s. of Eq. ( !T9|) with the replacement cos (20) Px = Pt cos 0. Figure [25] shows the re- 
sults for proton and pion sideflows in 10 AGeV Au+Au collision with 6 = 4 fm. On can see, 
that the calculation with the EOS-PT predicts a back bending in the central bin, \y\ < 0.5, 
which is not seen in data. On the other hand, the observed maxima of | < > | are stronger 
overestimated in the hadronic scenario. 

Figure [26] represents proton and charged pion elliptic flows for the same reaction as in 
Figs. [2H4251 In agreement with discussion in Sec. IIVB| the deconfinement phase transition 
leads to enhancement of elliptic flows as compared to the purely hadronic scenario. The 
relative increase of V2 is larger for pions. As before, we consider two values of the impact 
parameter, 6 = 4 and 7 fm. Experimental data exist only for They are better repro- 
duced with the EOS-HG. However, strong sensitivity of results to the choice of centrality 
and possible influence of dissipative effects (see below) do not permit drawing a more definite 
conclusion. 
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FIG. 26: (Color online) Elliptic flows of protons (left panel) and charged pions (right panel) as 
functions of rapidity in 10 AGeV Au+Au collisions. Solid (dashed) lines are calculated with EoS- 
PT (EoS-HG). Thick and thin lines correspond to 6 = 4 and 7 fm, respectively. Full dots are 



experimental data 



69( 1 for negative pions in 10.6 AGeV midcentral Au+Au collisions. Circles are 



obtained by reflection. 

To check the sensitivity of collective flows to the bombarding energy, we have performed 
calculations for Pb+Pb collisions at i^iab = 40 AGeV. The results are shown in Figs. [2714281 
As compared to the AGS energy (see Fig. the proton directed flow is not so sensitive to 
the EoS, at least in the central rapidity region. As one can see from Fig. [271 fhe calculations 
with the EoS-PT again predicts formation of the pion antiflow. In this case experimen- 
tal data are in qualitative agreement with our calculation while the EoS-HG leads to the 
opposite (positive) slope of the directed flow. 

Similar conclusion about smaller sensitivity of collective flows to the phase transition 
at SPS energies may be drawn from Fig. [281 Here we show the rapidity dependence of 
proton and pion elliptic flows for 40 AGeV Pb+Pb collisions. Note that experimental values 
of V2 strongly depend on the methods used for their determination. One can see that the 
model overestimates the proton data at large cm. rapidities. As discussed in Ref. |55| . 
more definite conclusions can be made only after dissipative and freeze-out effects will be 
carefully evaluated. 
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FIG. 27: (Color online) Directed flows of protons (left panel) and charged pions (right panel) as 
functions of rapidity in 40 AGeV Pb+Pb collisions. Solid (dashed) lines are calculated with EoS- 
PT (EoS-HG). Thick and thin lines correspond to 6 = 4 and 5.6 fm, respectively. Full symbols are 
measured experimental data [7C|], while the open ones are obtained by reflection. 
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FIG. 28: (Color online) Same as Fig. [271 but for elliptic flows of protons (left panel) and charged 
pions (right panel). 
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VII. CONCLUSIONS AND OUTLOOK 



We have developed a (3+1) dimensional ideal hydrodynamical model to study heavy-ion 
collisions in the energy range 1-160 AGeV. The sensitivity of collective flow observables and 
particle spectra to the flrst-order deconfinement phase transition have been investigated. 
Our analysis shows that maximal values of energy- and baryon densities in the central box 
of the colliding system are signiflcantly larger if this phase transition occurs at some inter- 
mediate stage of a heavy-ion collision. As compared with the purely hadronic scenario, the 
calculations with the deconfinement phase transition predict broadening of proton rapidity 
distributions, enhancement of elliptic flows and appearance of directed antiflow in the cen- 
tral rapidity region. It it shown that the collective flow parameters are especially sensitive 
to EoS at E'lab — 10 AGeV. Our calculations with the EoS-PT predict that at such energies 
the system spends the longest time in the mixed phase. Similar conclusion that this energy 
domain is optimal for searching signatures of the deconflnement phase transition in nuclear 
collisions has been made earlier in Ref. 71 1. 



Our analysis does not allow us to decide which EoS, with or without the phase transition, 
gives better agreement with experimental data. One should bear in mind that our results 
are obtained by assuming a rather simple isochronous freeze-out and neglecting dissipative 
effects. To detect clear signatures of the deconflnement transition in nuclear collisions, new 
detailed measurements of collective flow observables, especially in the central rapidity region, 
would be needed in NICA and FAIR experiments. In the future we plan to generalize our 
model by introducing dissipative terms and a more realistic description of freeze-out process. 
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